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Abstract. - We study the mechanisms responsible for quantum diffusion in the quasiperiodic 
kicked rotor. We report experimental measurements of the diffusion constant on the atomic 
version of the system and develop a theoretical approach (based on the Floquet theorem) 
explaining the observations, especially the "sub-Fourier" character of the resonances observed 
in the vicinity of exact periodicity, i.e. the ability of the system to distinguish two neighboring 
driving frequencies in a time shorter than the inverse of the difference of the two frequencies. 



Quantum chaos is the study of quantum systems whose classical limit is chaotic. A major 
challenge of quantum chaos is to understand the mechanisms that make quantum chaos dif- 
ferent from classical chaos. An important difference between classical and quantum systems is 
the existence in the latter of interferences between various paths. At long times, a large num- 
ber of complicated trajectories interfere. One could expect the contributions of the various 
paths to have uncorrelated phases, so that the interference terms vanish in the average after 
some time, implying that quantum and classical transport should be identical. This simple 
expectation is however too naive, because phases of the various paths are actually correlated; 
this is for example the case for the kicked rotor. 

The quantum kicked rotor has been extensively studied experimentally in recent years 
[1-4]. In its atomic version it consists of a cloud of laser-cooled atoms exposed to short pulses 
of a far detuned, standing laser wave, corresponding to the Hamiltonian (for the external 
motion of the atoms) 



where 9 is the 2ir— periodic position of the rotor, p the conjugate momentum, T the period 
and K is proportional to the strength of the kicks. In the classical limit, this system is 
chaotic for K > 1 [5], and the motion is an ergodic diffusion in momentum space for K > 5. 
Dynamical localization (DL) is the suppression of such a diffusion in the quantum system by 
subtle quantum interference effects [6]. 
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Fig. 1 - (a) Averaged squared momentum (p 2 ) (p is in recoil momentum units) vs. number of double 
kicks. All curves show DL after 5-10 kicks. After this break-time, the dynamics is frozen for the 
periodic system (r=l, solid curve), but a residual quantum diffusion is observed in the quasi-periodic 
case, for r = 1 + 0.00111 and r = 1 - 0.00111 (dotted curves) and r = 1 ± 0.00222 (dashed curves), 
(b) The quantum diffusion constant (in units of the classical diffusion constant) versus r, from data 
shown in (a). It (almost) vanishes at r = 1 because of DL, and rapidly increases symmetrically on 
both sides, displaying a triangular cusp. 

DL is expected to disappear if the time periodicity is broken [7] , as experimentally observed 
in [3] . This can be done simply by adding a second series of kicks with a different period rT, 
giving 

v 2 K 

H(r, A) = y + — cos 9 ^ [S(t - nT) + 5(t - nrT - AT)] (2) 

n 

where A is the initial phase between the two kick sequences. If r is rational, the system is 
strictly time-periodic and DL takes place, but it is rapidly destroyed around any rational 
number. One way of characterizing this sensitivity to the time-periodicity is to measure the 
average kinetic energy (p 2 ^ of the system as a function of time (or number of kicks), which 
allows to deduce its diffusion constant. In this paper, we present new experimental results 
and a theoretical interpretation of the physical mechanism responsible for the destruction of 
DL in the vicinity of rational r. 

Our experimental setup is described in [8] . Cold cesium atoms are produced in a magneto- 
optical trap, the trap is turned off, and a series of short pulses of a far-detuned (13.5 GHz ~ 
2600 r) standing wave (around 65 mW in each direction) is applied. At the end of the pulse 
series, pulses of counter-propagating phase-coherent beams perform velocity-selective Raman 
stimulated transitions between the hyperfine ground state sublevels F g = 4 and F g = 3. 
A resonant probe beam is then used to measure the fraction of transfered atoms, which 
corresponds to the population in a velocity class. Repeated measurements allow to reconstruct 
the atomic momentum distribution P(p). In the periodic case, one observes, for t > ti, where ti 
is the so-called localization time, two manifestations of DL: i) P{p) is frozen in a characteristic 
exponential shape P(p) ~ exp(— \p\/t), where £ is the localization length, and ii) the average 
kinetic energy of the system tends to a constant value, or, equivalently, the diffusion constant 
D = lim^oo (p 2 ) jt vanishes. 

Experimental measurements of (p 2 ) and the corresponding diffusion constant are displayed 
in Fig. d In the periodic case r = 1 (solid line), after an initial linear increase, the kinetic 
energy saturates to a constant value - corresponding to D(r = 1) rj in Fig. db). Due 
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Fig. 2 - Sub-Fourier resonances. The measured atomic population in the zero-momentum class, 
after 20 (upper solid curve) and 100 (lower solid curve) double kicks is plotted vs. the ratio r of 
the periods. The period of the first sequence is T — 27.8/is, and K ~ 10. At r = 1, dynamical 
localization is responsible for the large number of zero-momentum atoms. Away from the exact 
resonance, localization is progressively destroyed. Note the narrowness (5 times smaller than Fourier 
limit) and the triangular shape of the resonance line. Dashed lines are fits to the experimental lines 
using Eq. Q. 

to residual spontaneous emission by the kicked atoms [9], D(r = 1) is not strictly zero. In 
the quasiperiodic cases, there is a residual diffusion due to the quasiperiodicity, and Fig. ^b) 
shows that D{r) oc \r — 1|. 

Another way to characterize the observed residual diffusion is to measure the population P 
of the zero-momentum class ( 1 ), as a function of r. Fig.[2]shows the experimentally measured 
Po after 20 and 100 double kicks. This "localization resonance" displays a sharp peak at r = 1, 
indicating the presence of DL, and decreases rapidly on both sides, evidencing the residual 
diffusion. The plot presents two surprising features: (i) The resonance is very narrow: after 
N kicks, it could be argued that the two quasi-periods can be distinguished only if they differ 
by 1/N (in relative value). This would predict a width of the order of Ar = 1/N = 0.01 
for 100 kicks, whereas we experimentally observe a width five times smaller, 0.0018 ( 2 ). (ii) 
The "sub-Fourier" resonance is not smooth, but has a marked cusp at the maximum. The 
understanding of the underlying mechanism behind the quantum behavior of the system shall 
also allow us to explain these features. 

Our analysis of the residual diffusion takes the periodic case as the reference system, 
because, being periodic, it can be analyzed using the Floquet theorem. A Floquet state (FS) 
\ifk) is defined as an eigenstate of the unitary evolution operator U(T) of Hq over one period 
T: U(T)\ipk) = exp(— iefc) \ipk) where eu is called the eigenphase. The temporal evolution 
of any state |^>) after n periods is \^{nT)) — J2k Ck e~ mek \(fk) with Cfc = (y>jt|^>(0)). The 
evolution of any quantity can be calculated using the basis of FS, for example: 



The FS of the chaotic kicked rotor are well-known: they are on the average exponentially 

( 1 ) As the total number of atoms is constant, measuring P(p = 0) is equivalent to measure the width of the 
momentum distribution: P(p = 0) oc (p 2 ) • 

( 2 ) Sub-Fourier factors of 1/37 have been reported in [10], and factors of 1/60 have been experimentally 
observed by us. 
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Fig. 3 - Eigenphases et of the evolution operator over one period for the doubly-kicked rotor, corre- 
sponding to the experimental conditions, versus the phase A. The finite duration of the kicks (800 ns) 
is taken into account. We have only plotted initially populated Floquet states having a significant 
weight: |{i/>(0)|v?fc(Ao)}| 2 > 10~ 4 (thin line), or > 10~ 2 (thick line). States appear (disappear) as 
their weights go above (below) the threshold. Note that Floquet states rapidly change when avoided 
crossings with other states are encountered. 

localized in momentum space around a most probable momentum p^, with a characteristic 
localization length I [11, 12]. Such a localization - from which DL originates - is far from 
obvious and is closely related to the Anderson localization in time-independent disordered 
one-dimensional systems [13]. 

The initial state is supposed to be localized in momentum space around zero-momentum, 
with a width much smaller than the width I of the FS (which is the case in the experiment). 
Hence, only FS with roughly \pk\ < I will play a significant role in the dynamics; we shall call 
such states initially populated FS. Eq. © is a coherent sum over FS. However, as times goes 
on, non-diagonal interference terms accumulate larger and larger phases. In a typical chaotic 
system, these phases will be uncorrelated at long enough times, leaving an incoherent sum 
over FS: 



This equation is valid once DL is established (i.e., for t > ti). How long does it take for 
the phases n{e% — ey) to be scrambled? This can be simply estimated from the level spacing 
between initially populated FS, and turns out to be roughly tg — IT, while (p 2 ) saturates to 
a value oc £ 2 . 

A key point is to realize that, if r is close enough to 1, the sequence of kicks is very sim- 
ilar to a periodic doubly- kicked rotor for r = 1, except that the phase A between the two 
sequences slowly drifts along the sequence. A small part of the kick sequence around the 
n th kick will seem "instantaneously" periodic, with a phase difference A(n) = Ao + n(r — 1), 
where Ao is the initial phase between the two kick sequences. The evolution operator of 
the quasiperiodic system from time [n — 1)T to nT is thus given by the evolution operator 
17(A) = exp(-ip 2 AT/ft) exp(-iK cos 6/hT) exp(-ip 2 (l - \)T/h) e?q>(-iK cos 6/ hT) of the pe- 
riodic doubly-kicked rotor. The total evolution operator can thus be written as a product of 
the "instantaneous" evolution operators n«==i ^(^o + n ( r ~ •"•))• For smai l enough \r — 1|, the 
adiabatic approximation [14] applies: if the system is initially in a FS of U(Aq), it remains 




(4) 



H. LlGNIER et al: QUANTUM DIFFUSION IN THE QUASIPERIODIC KICKED ROTOR 



5 



in the corresponding FS of the "instantaneous" evolution operator as A changes. This is il- 
lustrated in Fig. |3 which shows the Floquet spectrum of the periodic doubly-kicked rotor, 
obtained by numerically diagonalizing U{\). When A is varied the eigenergies evolve along 
complicated "spaghetti", characteristic of quantum-chaotic systems, with a large number of 
avoided crossings (AC). 

In the experiment, the initial state is a linear combination of FS. The adiabatic approxi- 
mation implies that it remains a linear combination of the "instantaneous" FS with the same 
weights (the phases evolve, but the squared moduli remain constant). As discussed above, the 
coherences between FS vanish after fy, which implies that Eq. JIJ remains valid close to the 
resonance, provided one uses the "instantaneous" Floquet eigenbasis: 

<p 2 (A))«5> fe | 2 (^(AJIpVM). (5) 

The adiabatic approximation implies that there are no population exchanges among Floquet 
eigenstates. This means that the weights |cfc| 2 are constant all along the evolution; they are 
evaluated at the initial time, corresponding to the initial wavefunction |^(0)) and to A = Aq: 
|cfc| 2 = |(?/>(0)|</?fc(Ao)}| 2 , where \<fik(^o)) are the "instantaneous" eigenstates corresponding 
to the initial time. The evolution of the average kinetic energy (p 2 ) is thus entirely due to 
the evolution of the "instantaneous" Floquet spectrum with the parameter A, which evolves 
adiabatically from its initial value Ao to the final value A = A + N(r — 1) corresponding to 
the end of the kick sequence. Thus, two distinct types of localization properties come into 
Eq. (JSJ): those of the initial Floquet spectrum, present in the constants \ck\ 2 , and those of the 
"instantaneous" Floquet spectra, present in the "instantaneous" eigenstates \ipk(X)). The |cfe| 2 
coefficient gives important weights to the FS represented in the initial distribution (that we 
shall call "initially populated" FS), which, because the initial momentum distribution is very 
sharp compared to the localization length, are mostly localized around zero momentum. This 
limits the range of the sum in Eq. J5j to eigenstates |y?fc(Ao)) centered at momenta |pfe(Ao)| < I. 
As A moves away from Ao, the center p^(A) of \(pk(\)) moves away from zero momentum. As 
these eigenstates keep the same weight in the sum, this enlarges the momentum distribution. 
This is the mechanism at the origin of the quantum diffusion responsible for the growth of 
the average kinetic energy in the quasi-periodic case. One thus expects (p 2 (X)) to have a 
minimum at A = Ao and to rapidly increase as the kicks are applied, that is, as A evolves. 
Since A — Ao is proportional to r — 1, (p 2 ) shall also present a sharp minimum at r = 1, which 
implies that the population in the zero-momentum class shall present a sharp maximum at 
r = 1, as experimentally observed in Fig.0 

What determines the linewidth of the resonance? A FS may considerably change because 
of its AC with other FS, as tiny AC may be crossed diabatically. On the average, the typical 
variation A — Ao of the phase between kick series for which \ipk(X)) loses the localization 
property of |<£>fe(Ao)}, is the distance AA C to the next AC. We immediately deduce that the 
full width A/ (in frequency) of the sub- Fourier line is such that: 

2AA C = (NT)Af. (6) 

If the classical dynamics is regular, the Floquet eigenstates evolve smoothly with the parameter 
A; a change in A of the order of one is thus required to significantly modify the Floquet states: 
2AA C w 1, which corresponds to A/ = I /(NT), that is, the Fourier limit. In a classically 
chaotic system, however, the level dynamics displays plenty of AC, see Fig. 01 AA C is then much 
smaller than unity, leading to sub- Fourier resonances. Eq. @ also predicts the linewidth to 
be inversely proportional to the temporal length of the kick sequence beyond the localization 
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time (i.e. the sub- Fourier character is independent of TV), as numerically observed [10]. The 
critical value AA C depends on the detailed dynamics of the system. It can be roughly estimated 
by visual inspection of the quasi-energy level dynamics, Fig- El to be of the order of 0.05 and 
a "factor 10" sub-Fourier line, about twice the experimentally observed factor (the additional 
experimental broadening is due to the transverse profile of the laser mode leading to spatial 
inhomogeneities in K). 

From Eq. JSJl and the previous analysis, it is expected that (p 2 ) significantly increases 
from its minimum value at A — Ao- AC between FS localized around the same momentum 
are rather large (this is what determines AA C ) whereas AC between states localized a distance 
L 3> t apart in momentum space are typically much smaller and scale like exp(— L/t), as a 
consequence of the exponential localization of the FS. There is thus a very broad distribution of 
AC widths, with a large number of tiny AC. A tiny AC typically extends over a small A interval 
and thus tends to produce small values of AA C . The increase of (p 2 ) with A thus depends on 
the number of small AC encountered. In the presence of exponential localization, the AC 
density scales with size C as 1/C for C — ► 0, and (p 2 (\)) — (p 2 (Ao)) shall behave like |A — Ao| 
( 3 ), producing the cusp experimentally observed in the resonance line, Fig. The large 
number of extremely small AC is responsible for the singularity of the sub-Fourier resonance 
line. Another consequence is the diffusive behavior observed in the vicinity of the resonance, 
see Fig. nfb). Indeed, (p 2 (A)) — (p 2 (Ao)) increases linearly with |A — Ao|, itself a linear function 
of time and of \r — 1|. Thus, our model correctly predicts two non-trivial properties: {p 2 ) 
increases linearly with time and the corresponding diffusion constant is proportional to \r— 1|. 
This is distinct from the prediction of Random Matrix Theory, which implies an increase of 
(p 2 ) initially quadratic in A — Ao [15]. 

There will always be some degree of nonadiabaticity. Whatever small r — 1| is, tiny enough 
AC will be crossed diabatically. This puts a lower bound on the size of the AC effectively 
participating in the quantum dynamics and produces a rounding of the top of the sub- Fourier 
resonance line, too small to be seen in the experiment after 100 kicks, but easily visible after 
20 kicks, Fig. □ 

Our approach concentrates on the immediate vicinity of the resonance. What happens 
in the wings of the sub-Fourier line? Eq. JSJ indicates that this depends on the residual 
correlation between ^.(Ao)) and \ipk(\)) for |A — Ao| > AA C . A quantitative answer to this 
question is difficult. However, it seems clear that the quantum diffusion constant does not 
exceed the classical one, which corresponds to vanishingly small interference terms. Random 
Matrix Theory tells us that this type of parametric correlation usually decays algebraically 
with A — Ao- We thus propose the following ansatz: 

(p>T)) = (p^ + ^ |r _^^ AA/ T (7) 

where D c i is the classical diffusion constant and (p 2 )dl the saturation value of p 2 due to DL. 
This equation fits (using P(p = 0) oc (p 2 )~ 1//2 ) very well the experimental curves in Fig. [21 
reproducing the linear behavior at the center of the resonance and the classical diffusion in 
the wings. The parameters are however fitted, not extracted from the previous analysis. The 
reason for that is that, in the experiment, the standing wave laser beams creating the kicking 
potential have a gaussian profile whose width is only a few times the width of the spatial 
atomic distribution. This means that atoms away from the center of the beam see a smaller 
value of K than the atoms at the center of the beam. There is thus an averaging effect over 



( 3 )This basically comes from the fact that the probability to encounter a "bad" AC which expels the atoms 
from the zero- momentum region is proportional to the length of the interval |A — Ao |- 
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K to which the dynamics is sensitive, unfortunately preventing us form directly comparing 
the fitted values to the theoretical values of the parameters. 

In summary we have developed a theoretical approach for the mechanism of sub-Fourier 
resonances, which correctly predicts the unexpected observed features. In particular, our 
approach evidences the role of chaotic dynamics for producing narrow sub- Fourier resonances. 
Deviations from exact periodicity are treated in the framework of the adiabatic approximation 
for the Floquet spectrum of the system, which thus goes beyond a pertubative approach. The 
dynamics is governed by instantaneous Floquet eigenstates, which are non-trivial objects, as 
they need to have very well defined internal (i.e. between the various parts of the wavefunction 
for each Floquet state) phase coherence to be stationary states of the periodic system, but at 
the same time, the inter-state coherences do not play any role (beyond the localization time, 
Floquet states effectively add incoherently). The dynamics is dominated by an incoherent 
sum of internally extremely coherent states. It shows that the role of interferences in quantum 
mechanics is far from obvious, and can produce unexpected behaviors. 

Laboratoire de Physique des Lasers, Atomes et Molecules (PhLAM) is Unite Mixte de 
Recherche UMR 8523 du CNRS et de l'Universite des Sciences et Technologies de Lille. Lab- 
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